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Abstract: We consider dark sector particles (DSPs) that obtain sizeable interactions with 
Standard Model fermions from a new mediator. While these particles can avoid observation 
in direct detection experiments, they are strongly constrained by LHC measurements. We 
demonstrate that there is an important complementarity between searches for DSP produc¬ 
tion and searches for the mediator itself, in particular bounds on (broad) dijet resonances. 
This observation is crucial not only in the case where the DSP is all of the dark matter but 
whenever — precisely due to its sizeable interactions with the visible sector — the DSP 
annihilates away so efficiently that it only forms a dark matter subcomponent. To highlight 
the different roles of DSP direct detection and LHC monojet and dijet searches, as well as 
perturbativity constraints, we first analyse the exemplary case of an axial-vector mediator 
and then generalise our results. We find important implications for the interpretation of 
LHC dark matter searches in terms of simplified models. 
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1 Introduction 

In the event that the Large Hadron Collider (LHC) discovers evidence for new invisible 
particles from telltale missing energy signatures, it would not be possible to unequivocally 
establish if such a dark sector particle (DSP) is a viable dark matter (DM) candidate. To be 
able to draw such a conclusion, it is essential to combine LHC searches with cosmological 
or astrophysical observations as well as with results from direct and indirect detection 
experiments. Doing so, however, necessarily requires some assumptions on the properties 
and interactions of the DSP. 
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By interpreting potential missing energy signatures within the context of a specific 
theoretical model, it is possible to make concrete comparisons between collider searches, 
cosmology, and direct and indirect detection signatures [1,2]. A well-motivated and simple 
class of models involves an additional scalar or vector mediator in such a way that the 
DSP interacts with SM states via the exchange of this new particle. An attractive feature 
of this framework is that, in the parameter regions currently probed by experiments, the 
interactions can be large enough to allow for sizeable DM annihilation rates in the early 
Universe and hence to avoid DM overproduction [3]. One possible example model is a 
massive Z' arising from a new broken U(l)x gauge symmetry [4, 5], which can in principle 
have large couplings to both the DSP and SM fermions [6-12]. 

If this new s-channel mediator has a mass comparable to LHC energies, it will affect 
missing energy signatures at colliders and must therefore be taken into account for the 
interpretation of DM searches at the LHC and the comparison with direct detection ex¬ 
periments [8, 13-17]. Furthermore, the presence of such a mediator has two additional 
important implications, which motivate the central themes of this work: 

• As well as the usual searches for DM at the LHC based on missing energy in associa¬ 
tion with SM particles one also needs to consider dedicated searches for the mediator 
particles themselves, which make use of the fact that any mediator produced from 
SM particles in the initial state can also decay back into SM states. To illustrate the 
complementarity of these two approaches we will consider monojet searches [18, 19] 
as a typical representative of the former category and compare them to constraints 
on the mediator arising from dijet resonance searches. We will show that, by com¬ 
bining data from UA2 [20], the Tevatron [21] and the LHC [22-24], searches for dijet 
resonances probe a wide range of mediator masses from the electroweak scale to well 
above the TeV scale. Crucially, these searches turn out to be sensitive even for broad 
resonances, thus ruling out large swathes of parameter space which would otherwise 
appear viable when confronted with only monojet and direct detection constraints. 

• In the parameter regions probed by LHC monojet and dijet searches, the generic pre¬ 
diction is that the DSP would be underproduced in the early Universe. In particular, 
if the DSP mass m x is larger than the mediator mass Mr, the process XX RR 
can easily deplete the DSP abundance, and for m x ~ Mr/2 annihilation into SM 
fermions will receive a resonant enhancement. This prediction is still consistent with 
all observations if it is assumed that the particle under consideration only constitutes 
a DM subcomponent and another particle (for example the axion [25, 26]) makes up 
for the remaining DM abundance. 1 As a result, the local density of DSPs in the Milky 
Way can be significantly smaller than what is usually assumed for the interpretation 
of direct detection experiments. The predicted DSP underabundance must therefore 
be taken into account when comparing between direct detection and LHC constraints 

1 An alternative approach would be to call the standard thermal history of the Universe into question 
and consider additional dilution or production mechanisms for DM [27]. Such modifications would however 
require further fields and couplings. In the present paper, we therefore take the standard thermal history 
as a fundamental underlying assumption. 
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Figure 1. The processes considered in this work in terms of visible sector quarks (q, q ), DSPs (%, \) 
and the on-shell (off-shell) mediator particle R (R*). The various process are: (a) DM annihilation 
which sets the relic abundance, (b) DM scattering in direct detection experiments, (c) monojet 
signatures, in this case due to initial state radiation of a gluon, (d) LHC Dijet resonance signatures 
purely through mediator-quark couplings and (e) dijet associated production. 


in order to avoid overstating the strength of direct detection limits. This approach 
leads to a compelling interplay between the different DM detection techniques and 
will lead us to conclude that the LHC monojets, LHC dijets and direct detection 
strategies each has a unique foothold in the search for DSPs. 

In figure 1 we sketch the setup for a dark sector theory involving a DSP x and a 
mediator between the visible sector and the dark sector R, together with the detection 
processes considered in this work. We denote the couplings between the mediator and 
the visible sector quarks (the DSP) with g q (g x ). For the purposes of exploring the broad 
phenomenology of this dark sector and the general interplay between the different probes let 
us combine the two couplings into an effective DSP-SM coupling g = <Jg q g x and consider 
the effect of varying the coupling g. The local density of DSPs in the Milky Way p is 
proportional to the DSP relic abundance from thermal freeze-out Hdsp, which scales as 
the inverse of the annihilation cross section, i.e. p oc Hdsp oc g~ 4 . Any cross section 
involving interactions between the visible sector and the DSP, such as collider production 
and direct detection, will scale as a oc g 4 [1, 28-31] (assuming an off-shell mediator). Thus, 
broadly speaking, the rate of events in different DM probes have very different scaling with 
couplings if a standard thermal history is assumed. They are: 

• Collider searches for missing energy: Rate oc a oc g 4 . 

• Direct detection: Rate oc (a x p) oc g° . 

• Indirect detection: Rate oc (a X p 2 ) oc g ~ 4 . 

Furthermore, resonance searches at colliders typically depend on the production cross sec¬ 
tion for the resonance, or, multiplied with the branching ratio into the final state under 
consideration. If the (on-shell) mediator has a large branching into light quarks we hence 
obtain the final important signature 

• Collider searches for dijet resonances: Rate oc ctr oc g 2 . 

This simple consideration demonstrates that, assuming a standard thermal history and con¬ 
sidering the specific phenomenology of the mediator, these four different detection strate¬ 
gies are parametrically complementary. In essence, large couplings imply large collider 
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rates but small abundances, small couplings imply small collider rates and large abun¬ 
dances. Interestingly, to a first approximation, direct detection of a DSP is independent 
of the interaction strength g. Another pertinent consequence is that the dijet constraints 
break the degeneracy between mediator couplings to the visible and dark sectors. 

Of course, the discussion above has been intentionally over-simplified. In reality even 
in the simplest scenarios there are four relevant parameters, namely m x , Mr, g x and g q . 
Specifically, detection strategies often exhibit a highly non-trivial dependence on the medi¬ 
ator mass Mr and also on the ratio m x /MR, which determines whether the mediator may 
decay to DM, whether there can be resonant enhancements and whether the annihilation 
channel XX Ah’ is open. In this work we will take all of these parameters into account 
and furthermore consider different coupling structures for g q in order to thoroughly explore 
the relevant phenomenology and interplay between different dark sector probes. 

Moreover, a model containing just a DSP and a mediator is typically only a low-energy 
effective description of a more complete theory and we therefore have to pay attention to 
the validity of our approach. In particular, we have to ensure that the model remains 
perturbative and does not violate unitarity at the energies at which it is being probed [32- 
36]. As we will show, these considerations imply that, for the case of an axial-vector 
mediator, the DM mass must not significantly exceed the mediator mass. The resulting 
theoretical constraints may not be as concrete as experimental exclusion limits, but must be 
kept in mind when studying the cosmology and collider phenomenology of a DSP interacting 
via a mediator. 

In the present work we will mostly focus on the scenario of an axial-vector mediator for 
concreteness and also comment on the case of vector couplings. This choice is motivated by 
the desire to demonstrate the complementarity between different approaches as clearly as 
possible as no single probe dominates the experimental limits for this model. Specifically, we 
will demonstrate that only small, often closed, parameter spaces survive current constraints, 
providing clear targets for future experimental exploration. Many of our conclusions apply 
also to other models of thermal DM that contain a new mediator. The examples chosen 
here are therefore meant to illustrate how this entire class of dark sector theories can be 
probed by searches across the frontier of dark sector physics. 

The paper is structured as follows. In section 2 we introduce the general interactions 
of the DSP and the mediator and discuss the coupling structures and parameter ranges 
relevant for our study. Section 3 then focusses on the calculation of the DSP relic abundance 
and the resulting rescaling of direct detection bounds. Monojet and dijet constraints will 
be discussed in section 4 and section 5, respectively. We present the combination of all 
of these constraints in section 6 and a discussion of more general scenarios and future 
prospects in section 7. 

2 Simplified interactions of a vector mediator 

We consider a simplified model of a DSP %, taken to be a Dirac fermion, and a vector 
mediator R. In this model no gauge symmetry, additional fields to cancel anomalies or 
symmetry-breaking structure is specified. Thus, as with any effective theory, new dynamics 
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will enter to complete this model in the ultraviolet. We assume that at low scales the 
specific ultraviolet-completion is decoupled so that x and R are the only fields relevant for 
our study. The couplings of R and x are written as 

£ds = xY'(9x - 9$ 7®)X • (2-1) 

In the same way, the couplings of R to SM particles are given by 

£ // = E R * f -- 9f 7 5 )/ , ( 2 . 2 ) 

f=q/ 

where q and t denote SM quarks and leptons, respectively. The partial decay widths of the 
mediator in terms of these couplings are 

r (R -> XX) = )() y 4 [(sj) * 2 + (si 1 ) 2 + z x ( 2 (sf) 2 - 4 (Sx ) 2 )] ’ 

nR^Sf)= i ^Xi^[{3j? + (gff + z,(2(a)f-i( g jf)} , (2.3) 

where z x j = m 2 f/M )|, N c = 3 for quarks and N c = 1 for leptons. 

If both g x and g q are non-zero, the DSP interacts with nuclei via spin-independent 
(SI) scattering. These interactions receive a coherent enhancement proportional to the 
square of the target nucleus mass and are therefore strongly constrained by direct detection 
experiments [37] (see e.g. [38] for a recent discussion of these constraints in the context of 
vector mediators). In this case, LHC bounds typically give no relevant constraints unless 
the DSP is very light. 

For the largest part of this work we will therefore assume that at least one of the 
vector couplings vanishes. 2 In this case, DSP-nucleus scattering will be dominated by the 
spin-dependent (SD) interactions induced by the axial couplings. Potential cross-terms 
such as g q g x lead to momentum suppressed scattering in the non-relativistic limit, which 
can safely be neglected. 2 To simplify our analysis, we will set both g q and g x equal to 
zero unless explicitly stated otherwise. More general couplings structures will be discussed 
in section 6.3. 

We assume that there is no direct link between the quark couplings g q and the leptonic 
couplings gg, in contrast to the case in which the mediator R obtains its SM couplings 
from mixing with the Z boson. The reason is that the leptonic couplings are very tightly 
constrained by searches for dilepton resonances [6, 10, 11]. Rather than forcing the quark 
couplings to be equally small, we want to consider the case where the mediator couples 
much more strongly to quarks than to leptons (see e.g. [42, 43] for a discussion of a baryonic 
Z'). In this case, the leptonic couplings will not give a relevant contribution to the DM 
phenomenology of the model and we can simply set g\ = gf = 0. Moreover, it was shown 
in [8] that similarly stringent constraints apply to the interactions between R and SM gauge 
bosons as well as the Higgs boson, which are therefore also assumed to be negligible. 

2 For example, the vector coupling g^ automatically vanishes if the DSP is a Majorana fermion. 

2 Note that these interactions also induce SI but mass-suppressed scattering at one-loop level. The 
resulting contributions are, however, subdominant to tree-level SD interactions [39-41]. 
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Figure 2. Left: Total width of the mediator R divided by Mr as a function of the mediator 
mass for various coupling ratios g x /g^- Right: The corresponding branching ratio into dijets as a 
function of Mr. For both plots, the DSP mass has been fixed to m x = 100 GeV. 

In summary, we will focus on the following simplified interactions: 

X7 m 7 5 X R» + E 4 R„ ■ (2.4) 

9 

We allow for the possibility to have different couplings for the individual quark flavours. 
To study how the phenomenology of the model depends on the different couplings, we will 
consider a number of different scenarios: 

• Universal couplings: g£ = gj = gj = gj = g£ = gf = gf. 

• Isovector couplings: g£ = g£ = gf = -gf = -gf = -gf = gf. 

• No heavy-quark couplings: gf = gf = gf = gf = gf and gf = gf = 0. 

As we will see below, the heavy-quark couplings play a non-negligible role, because they 
can provide additional annihilation channels for the DSP and thereby reduce the total DSP 
abundance (thus weakening bounds from direct detection experiments) and because they 
can increase the width of the mediator (thus weakening bounds from LHC searches). 

Clearly, there is no fundamental reason why gf should be identical to gf. Even if both 
couplings arise from a new gauge group, the respective charges and hence the resulting 
couplings could easily differ by a factor of a few [44, 45]. Alternatively, if the DSP couples 
directly to the mediator and couplings to quarks are generated only through a small mixing, 
gf could easily be significantly larger than gf. In the following, we will usually consider 
the three cases gf = gf, gf = 4 gf and gf = 9 gf. 

Since for processes involving an off-shell mediator, cross sections are typically propor¬ 
tional to (gf) 2 (gf) 2 , it will be useful to consider the effective coupling g = (gf gf) 1 / 2 . In 
figure 2 we show the total width of the mediator, as well as the corresponding branching 
ratios into jets, as a function of Mr for various coupling ratios gf/gf. For Mr 3> m x we 
find that the mediator decays dominantly into quarks if gf /gf < 4 and dominantly into the 
DSP if gf/gf > 4. If the coupling ratio is approximately gf/gf ~ 4 both decay channels 
have comparable branching ratios, leading to the smallest mediator width for fixed g. 
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To conclude this section, let us briefly discuss the parameter ranges of interest. First 
of all, the magnitude of the couplings is bounded from above by the requirement of per- 
turbativity, which implies g^ x < s/~A/. In addition, we require that the mediator width 
satisfies Tr/Mr < 1, such that it makes sense to interpret the mediator as a new particle. 
This second requirement implies roughly g// <1.5 and g x < 6. In practice, we will focus 
on the two cases g = 1 and g = 0.5, which for the coupling ratios discussed above give 
Yr/Mr < 0.5. We do not consider the case Yr/Mr > 0.5, because such mediators would 
typically be very difficult to observe in searches for dijet resonances. 

For the mediator mass, we restrict ourselves to the case Mr > 100 GeV. The reason 
is that we want to avoid the case where the mediator mass is very close to the mass of 
the Z-boson Mz, because for Mr ~ Mz quark loops induce large mixing between the 
two gauge bosons, leading to conflicts with electroweak precision observables [5, 46-48]. 
Parametrically, the kinetic mixing £ D e d ,J Z v (d^R v — d^R/j is given by 

' ~ E ~ 1 V) 2 < 2 - 4 5 ) 

<? 


and the correction to the p parameter, A p = M^/(M § cos 2 Q\ y) — 1, is roughly 


A p ~ e 2 


Mf 

M 2 - M| ’ 


( 2 . 6 ) 


which is sufficiently small for g^ < 1 and Mr > 100 GeV. 1 

Finally, if the theory is chiral, i.e. if g x / 0, the DSP mass cannot be raised arbitrarily 
compared to the mediator mass. The reason is that in a chiral gauge theory we essentially 
have to invoke a Higgs mechanism to generate both the DSP mass and the mediator mass: 


£ A — (D /i 4>) t - (y x 4> XLXR + h.c.) , 


(2.7) 


where 4* is the Higgs field in the dark sector and D^ = <9^ — ig'q R^ with g' and q being the 
gauge coupling of the new U( 1) / and the corresponding charge. In order for the Yukawa 
interaction to be gauge invariant, we must require q$> = q XL — q XR . Once 4> aquires a 
vacuum expectation value, (4>) = v'/V%, we obtain the mass terms 

£ A - ^g 2 q% v' 2 R^Rf, - ~^=y x v %X , (2.8) 

implying Mr = g' q<$, v' = g' (q XL - q XR ) v' = g£ v' and m x = y x v'/s/2 = y x Mr/ (s/ 2 g^). 
In order for the Yukawa interaction to remain perturbative we impose the bound 

m x < ^pr M R ■ ( 2 - 9 ) 

4 Another concern are loop-induced couplings of the mediator to leptons via the Z-mixing, leading to 

potential signals in searches for dilepton resonances. We have checked that these signals give weaker 

constraints than searches for dijet resonances. 
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We note that, even if the mediator mass and the DSP mass are generated in a different 
way, a similar constraint must apply to ensure that the coupling between the longitudinal 
component of R and the DSP remains perturbative. Indeed, we will see below that this 
inequality is closely related to the requirement of perturbative unitarity in DSP scattering 
and annihilation. 


3 Non-collider constraints 
3.1 Relic density 

For the Lagrangian introduced above, the DSP annihilation cross section into quarks is 


™{xx -t qq) - 


3 m 2 yj 1 - z q 


2tt [(M 2 ~Amlf + {T R M R y 


(. 9x ) 2 [(9q ) 2 ( 2 + Zq) + 2 (gf ) 2 (1 - Zq)\ 




(4m 2 - MV) 2 


M 4 r 


(3.1) 


where z q = m 2 /m 2 and we have chosen unitary gauge for the propagator of the vector 
boson as this should capture the effects of the Goldstone bosons present in an ultraviolet- 
complete realisation of this model. 

If g[ is zero, the annihilation cross section is proportional to m 2 /m 2 , i.e. there is 
a helicity suppression for annihilation into light quarks. In this case, it is important to 
also include the p -wave contribution for calculating the DSP relic abundance. Expanding 
av = a + bv 2 , we find (setting g^ = g x = 0) 


b 


(<4) 2 «) 


2 m 2 x 


2vr [(Af£ - 4m|)2 + (r* M R f 


-(1 - z / f ' 2 • 


(3.2) 


For m x > Mr the DSP can also annihilate directly into the mediator, which then 
subsequently decays into SM particles. The corresponding annihilation cross section is 
given by 


(m 2 - Mj,) 3 / 2 

- RR) = 4 ^ 404 - 2 ™^ 

x (8(9^) 2 (sf ) 2 m 2 x + [(9^) 4 - HSxf(Sx) 2 + Of) 4 ] Mr) ■ (3-3) 

To calculate the relic density of the DSP, we have implemented both types of in¬ 
teractions in micrOMEGAs v3 [49], which numerically solves the Boltzmann equation. As 
micrOMEGAs takes into account the full expressions for the annihilation cross section (rather 
than the velocity expansions shown above), the calculation gives reliable results also close to 
resonance. We show the results of this calculation in figure 3. The left panel shows contours 
of constant relic density for fixed couplings g x = g' 4 = 1. The right panel shows the param¬ 
eter region excluded by the requirement Ddsp < Udm ~ 0.119 /h 2 for different choices of 
the coupling product g = (g£ g^) 1 / 2 , keeping the coupling ratio fixed to g x /g q = 1. In the 
approximation that Ddsp is proportional to g ~ 4 this plot contains the same information as 
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Figure 3. The relic abundance ^dsp of the DSP as a function of m x and Mr. Left panel: Contours 
of constant SIqsp f° r 9 X = 9 q = 1- Right panel: The parameters giving Odsp = Ddm for different 
choices of couplings. The two plots are not equivalent since the mediator width and the freeze-out 
temperature Tf depend on the couplings. 

the one on the right-hand side. For example, the line Odsp = 1CP 1 J)dm f° r 9 = 1 in the 
left panel corresponds to the exclusion line ^dsp = ^dm for g A = 0.1 in the right panel. 
Visible differences arise, however, from the fact that Odsp depends on the mediator width, 
as well as on the freeze-out temperature, which can be significantly different from the usual 
choice Tf ~ m x /20 in the case that Odsp is very different from Odm- 

We would like to emphasise that — strictly speaking — the calculation performed 
above provides only an upper limit on the relic abundance of the DSP, which could be 
further reduced in the presence of additional annihilation channels. As argued in section 2, 
however, experimental constraints essentially require the DSP to couple much more weakly 
to leptons and to SM bosons than to quarks, so that it appears difficult to obtain sizeable 
contributions from any of these interactions. In any case, such additional annihilation 
channels would only reduce the abundance of the DSP and therefore further weaken direct 
detection and indirect detection compared to LHC searches, so that we essentially consider 
the case that is most optimistic for direct detection.'’ 

Before we turn to the discussion of direct detection constraints, let us briefly return to 
eq. (3.1). Unitarity arguments require that the cross section behaves as av oc 1/s at high 
energies, which in the case of DM annihilation implies av oc 1 /m^. Thus in eq. (3.1) the 
term proportional to the axial couplings appears to violate unitarity in the limit m x —>• oo 
as the cross section scaling becomes av oc rn^/M^,. As discussed above, however, this 
limit is not physical, due to the diverging couplings between the DSP and the longitudinal 

5 It is rather difficult in general to raise the DSP abundance above the value predicted by naive thermal 
freeze-out. One attractive possibility would be to assume an initial particle-antiparticle asymmetry in the 
dark sector. In this case, the DSP relic abundance is essentially set by this asymmetry and does therefore 
not directly reflect the interaction strength of the DSP. We will return to this possibility in section 7. 
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component of R. Imposing the perturbativity constraint from eq. (2.9), rearranged to 
Mr > g^m x /VM, is sufficient to make sure that the annihilation cross section remains 
well-behaved with av oc 1 /m x for large DM mass, thus we will impose this constraint 
on the parameter space throughout. It would be interesting to investigate perturbative 
unitarity in this model further to see if this constraint could be strengthened. 


3.2 Direct detection 


As noted above direct detection experiments constrain the rate of interactions between the 
DSP and nuclei, which is proportional to the local DM density p. Using the results shown 
in figure 3, we can therefore now rescale the bounds from direct detection experiments 
under the assumption that p is proportional to Ddsp- (> To obtain constraints we use the 
recent results from the LUX experiment [37]. 

The SI and SD DSP-nucleon scattering cross sections at zero momentum transfer are 
given by 


si _ 

a N ~ JN 


Vnx 

7T M r 


(3.4) 


and 


,SD _ „2 3 ^Nx 
7 tM R 


a N ~ a N 


(3.5) 


where N stands for either p or n and p,Nx = m x rriN/(jn x + m n) is the reduced DSP- 
nucleon mass. The effective DSP-nucleon couplings for SI interactions are given by f p = 
g^(2g^ + g^) and f n = g^(gu + 2gJ). For SD scattering, we can write 


a P ,n = fJ x A '? (p,n) 9q ■ 

q=u,d,s 


(3.6) 


The coefficients Aencode the contributions of the light quarks to the nucleon spin, 
which can be extracted from polarised deep inelastic scattering. The Particle Data Group 
values are [50] 

Au (p) = A d {n) = 0.84 ± 0.02 , 

Ad (p) = Au (n) = -0.43 ± 0.02 , (3.7) 

As (p) = As (n) = -0.09 ± 0.02 . 


For most of the relevant parameter space the strongest constraints on SI and SD 
interactions come from the recent LUX experiment [37]. However, a dedicated LUX analysis 
for SD constraints has not yet been performed and thus we calculate our own constraints, 
similarly to [12, 14, 51]. We use the spin structure functions from [52] and the same limit¬ 
setting procedure as described in [53] to determine the parameter region excluded at 95% 

f ’This is a good approximation in the case that all of the DM is cold and collisionless. If, on the 
other hand, the dark sector contains a mix of warm and cold DM or contributions from DSPs with large 
self-interactions or dissipation, the scaling can become more complicated. We do not consider this case 
further. 
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Figure 4. Parameter regions excluded by LUX. Left panel: The conventional LUX bound for 
Udsp fixed to J7 dm for = g£ = 1 (purple) compared to the relic density constraint (red) and 
the rescaled LUX bound (orange). The perturbativity bound from eq. (2.9) is shown in grey. Right 
panel: The rescaled LUX bound for different coupling choices. 

confidence level (C.L.). This procedure leads to good agreement with the published LUX 
bound for the case of SI interactions in the DM mass range that we consider. 

The resulting constraints are shown in figure 4. The left panel illustrates how the 
direct detection bound is constructed. The grey region indicates the parameter space 
where the perturbativity bound from eq. (2.9) is violated. The purple shaded region shows 
the parameter space conventionally excluded by LUX when assuming Q = Ddm for g^ = 
g£ = 1. The red shaded region shows the parameter region excluded by the requirement not 
to overproduce the DSP. Finally, the orange region shows the parameter region excluded 
by the rescaled direct detection bound (allowing both Ddsp < Ddm and Ddsp > ^dm) 
obtained by taking the thermal relic density as a prediction of the theory. Clearly, in the 
red shaded region (where Ddsp is larger than Ddm) the rescaled bound is stronger than the 
conventional bound, while outside of this region it is weaker because the DSP constitutes 
only a fraction of the DM abundance. By construction, both bounds coincide along the 
line Ddsp = ^dm- 

The right panel of figure 4 shows the rescaled direct detection bounds for a variety of 
different couplings. As expected, we find that (for fixed ratio of g£ and g^) the bounds 
are largely invariant under an overall rescaling of the couplings, because the change in 
the scattering cross section compensates the change in the DSP density. Nevertheless, 
there are small observable differences due to non-negligible changes in the mediator width 
and the freeze-out temperature. The largest changes are observed in the parameter region 
m x > Mr, mt, where the annihilation channels XX ~> tt and XX —> RR open up and 
introduce a dependence of the relic density on g^ t and the coupling ratio g^/g^. In 
summary, we see that by taking the reasonable assumption that the thermal abundance is 
a prediction of the model the constraints from direct detection on a particular theory are 
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substantially altered, and much larger regions of parameter space are allowed than if the 
assumption of 17 = 17 dm is imposed. 

3.3 Indirect detection 

Finally, we would like to point out that in the parameter region allowed by the requirement 
17dsp < 17dm there are typically no observable signals from DM indirect detection in 
any present experiment and in particular our model does not provide an explanation for 
the diffuse GeV-energy excess of gamma-ray emission from the Galactic Centre observed 
with the Fermi-LAT instrument [54-58]. The reason is that for rn x < rrit the s-wave 
contribution to the DM annihilation cross section is helicity suppressed, so that freeze-out 
is dominated by the p- wave contribution, which becomes unobservably small in the present 
Universe. As a consequence, the parameter region in which the DSP corresponds to all of 
DM corresponds to av <C 2.5 • 10 _26 cm 3 /s in the present Universe. For m x > mt, on the 
other hand, there are presently no indirect detection experiments sensitive to the thermal 
cross section av ~ 2.5 • 10 _26 cm 3 /s for 17dsp = 17dm- 

Of course, it is possible in our model to obtain significantly larger annihilation cross 
sections, but only at the expense of depleting the DSP abundance. Since annihilation 
signals depend on the square of the DSP density, making the annihilation cross section 
larger will reduce rather than enhance the magnitude of any indirect detection signal. 
Consequently, one can never exceed the signal strength expected in the case that the DSP 
constitutes all of the DM.' 

4 Monojet searches 

LHC searches for jets in association with missing transverse energy ($t) place strong 
constraints on the interactions between quarks and the DSP [18, 19]. These constraints are 
most easily interpreted in terms of contact interactions [2, 28, 29, 31, 33, 35, 36, 59-65], 
but have also been interpreted in terms of the exchange of a vector mediator [8, 13-17, 66]. 
Here we follow closely the analysis presented in [14] and find good agreement with the 
results presented there. 

If the mediator is forced to be off-shell in the production of the DSP (either because it 
is too heavy to be produced on-shell at the LHC or because Mr < 2 m x so that decays into 
the DSP are kinematically forbidden), the monojet cross section at the LHC depends on the 
couplings of the mediator according to a(pp —> jxx) oc g 4 . If the mediator can be produced 
on-shell (and provided that we can treat the resonance in the narrow-width approximation) 
the monojet cross section will be proportional to a(pp — > jxx) c* (<?y 4 ) 2 x BR(i? -» xx)- 
As long as we keep the coupling ratio g^/gf fixed, the branching ratios of the mediator 
are independent of g , leading to a(pp -» jxx) oc g 2 ■ 

In many realistic cases, however, the width of the mediator may become so large that 
the narrow-width approximation is no longer valid. Furthermore, for m x close to Mr/2 
there can be relevant contributions from both on-shell and off-shell mediators. In practice, 

Annihilation of DM subcomponents can however provide a plausible explanation for indirect detection 
signatures corresponding to annihilation cross sections below the thermal one. 
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Figure 5. Exclusion bounds from the CMS monojet search in the m x -Mn parameter plane for 
different coupling configurations. 

monojet cross sections can therefore depend on all relevant parameters (m x , Mr, g x and 
g^) in a non-trivial way. 

To derive exclusion limits, we consider the most recent monojet search from CMS at 8 
TeV, based on an integrated luminosity of 19.7fb _ [18]. This CMS analysis considers events 
with large amounts of provided there is a primary jet (j \) with transverse momentum 
Pt > 110 GeV and pseudorapidity \r]\ < 2.4. A secondary jet (j 2 ) with px > 30GeV is also 
permitted if the two jets are not back-to-back: \Acf>(j ], . 72 )! < 2.5. Events with high-p-p 
tertiary jets, electrons or muons are vetoed. The analysis considers a number of different 
requirements for the amount of in the range of 250-550 GeV. 

We simulate monojet events using the implementation of an axial-vector mediator in 
the Powheg-Box v2 [67]. We employ the MSTW2008 parton distribution functions (PDFs) [ 68 ] 
and set the renormalisation and factorisation scale p dynamically, choosing p = Hr /2 with 
Hr = \J m \x + Pt j 1 + Pt,] X ■ For showering and hadronisation we use Pythia v 6 [69]. It 
was shown in [67] that — in the presence of a veto on tertiary jets — including next-to- 
leading order corrections in combination with parton showering only leads to a modest 
change in the monojet cross section of less than 10%. We neglect this small enhancement 
here and generate all events at leading order. 

Across the entire parameter range that we consider, the strongest bound is obtained 
for the requirement ft? > 450 GeV . 8 For this particular cut, the CMS results exclude new 
contributions to the fiducial monojet production cross section in excess of 7.8 fb at 95% 
C.L. 

The resulting bounds for the different scenarios introduced in section 2 are shown in 

s For very small values of m x and Mr one would expect that a significantly weaker cut on $ T gives 
the strongest bound, due to the missing transverse energy spectrum being rather soft. However, due to an 
upward fluctuation in the data at low $ T and a downward fluctuation at high $, T the observed bound is 
stronger for tighter cuts. 
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figure 5. We find these results to be in full agreement with the expectations discussed above. 
In particular, we observe that for small mediator masses, the bounds become independent of 
the mediator width and hence depend only on the product g^. For small DSP mass, the 
bounds depend sensitively on the mediator width (this is most obvious from comparing the 
green and the orange dotted line). However, changing both the width and the production 
cross section while keeping g£ fixed leads to negligible changes in the limit of small DSP 
masses (as can be seen from comparing the green and the red dashed line). 

5 Dijet searches 

Searches for new resonances in dijet final states have been carried out in several experi¬ 
ments. At present stringent bounds are provided by UA2 [20], CDF [21], CMS [22] and 
ATLAS [23, 24] analyses, some of which have already been applied in theoretical studies 
to constrain scenarios with an additional U{1)' gauge group (see e.g. [44, 45, 70, 71]). In 
the present section we recast these experimental analyses in the context of an axial-vector 
mediator R. As we will see, these constraints are complementary to the ones described 
above, being more stringent for large values of g^ and small values of g^ (or large DSP 
masses). 

Interestingly, no single experiment provides the strongest bound across the entire pa¬ 
rameter space under consideration. The reason is that for small mediator masses the QCD 
background, produced mainly by two gluons in the initial state, increases much faster with 
increasing centre-of-mass energy than the signal. Consequently, at the LHC the signal is 
overwhelmed by QCD events in this parameter region. The latest ATLAS and CMS dijet 
analyses therefore focus mostly on the region with dijet invariant mass rrijj > 1 TeV, while 
UA2 and the Tevatron still provide competitive bounds for smaller masses. The impressive 
performance of the LHC, however, allows to produce for the first time heavy dijet reso¬ 
nances in association with other SM particles, such as Z and W bosons, providing a key 
opportunity to suppress the QCD background. Indeed, we will show that a recent analysis 
of such events in the context of technicolour [24] can be recast to give strong constraints 
on mediator masses below 300 GeV. 

We implement the experimental searches by means of several public codes. To gener¬ 
ate dijet events at parton level at leading order we use MadGraph v5 [72] with the model 
file produced by Feynrules v2 [73] and production at leading order using the CTEQ 6L1 
PDFs [74], 9 Parton level events are first passed through Pythia v6 [69] to simulate ini¬ 
tial and final state radiation, fragmentation and hadronisation, and subsequently through 
Delphes v3 [76] for fast detector simulation. 10 We finally employ Fast Jet v3 [77] for jet 
reconstruction, and MadAnalysis v5 [78] for cuts and data analysis. Part of the statistical 
analysis is carried out by means of MCLimit [79]. Details on each analysis are described 

9 In order to recast the dijet analyses, both the SM background and the observed data are taken from 
the experimental publications, while the signal is computed at leading order. Since next-to-leading order 
corrections generally lead to an enhancement of the resonance production cross section [75], our results are 
expected to be conservative. 

10 Note that, while detector effects give a small effect in monojet searches [19], they play an important 
role for the dijet analysis, because they can change the shape of the resonance. 
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Figure 6. Left: Reconstruction of the dijet invariant mass distribution for a narrow (orange) and 
a broad (blue) resonance with mass Mr = 200 GeV at UA2. The two vertical lines enclose the 
region between (R±2a)rrijj. Right: The solid lines indicate the dijet total production cross section 
times efficiency for different values of the mediator couplings with Mr = 200 GeV. The parameter 
choices corresponding to the blue (orange) line in the left panel are represented by a blue star 
(orange diamond). The dashed lines are obtained by using a constant efficiency of e = 70%, which 
is valid only for narrow resonances. 


below, starting from the ones relevant for the smallest mediator masses. 


5.1 UA2 dijet analysis 

The UA2 dijet analysis [20] is based on a data sample of 10.9 pb _1 of pp collisions at a 
centre-of-mass energy of y/s = 630 GeV. The analysis essentially requires events with two 
leading jets in the central region of the detector, with cos 6 < 0.6, and no further jets with 
transverse energy Er > 30GeV. A final cut rejects events with a dijet invariant mass, rrijj, 
outside the interval (R ± 2cj)Mr, with R ~ 0.95 and a ~ 0.085 (for the precise values of R 
and a, and their dependence on Mr, see [20]). As shown by the UA2 collaboration, this 
last cut has an efficiency e ~ 70% for dijet events coming from a narrow-width resonance. 
The analysis is sensitive to the region 130 GeV < rrijj < 260 GeV. 

Implementing this analysis in our Monte Carlo simulations, we check the efficiency e 
and find excellent agreement with [20] whenever Tr <C Mr. However, as expected, the 
efficiency can be drastically smaller for broader mediators. This effect is displayed in the 
left panel of figure 6 where we plot the (normalised) rrijj distributions for dijet events 
coming from a resonance with Mr = 200 GeV and universal couplings to quarks. In the 
figure the orange histogram corresponds to a narrow resonance, with g^ = 1/6, g x = 3/2 
and m x = Mr so that the invisible decay channel of R is kinematically closed. In this case 
the width of the mediator is about Tr/Mr ~ 0.01 and hence the width of the reconstructed 
distribution is dominated by final-state radiation and detector effects, which give rise to a 
displacement of the peak and tail towards the low-mass region. Despite this displacement, 
around 70% of the events fall within the region (RE 2a)M r (corresponding to the interval 
within the vertical dashed lines) in agreement with the UA2 finding. The blue histogram 
instead depicts the scenario of a broad resonance, with g£ = 1/3, g x = 3 and rn x = 40GeV. 
In this case the larger width of T r/Mr ~ 0.1 (convoluted with the proton and antiproton 
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PDFs) makes the resonance more likely to be produced at lower invariant masses rrijj. 
Consequently, less than half of the signal events fall within the (R ± 2ct)Mr region, and 
the efficiency e decreases to about 40%. 

The dependence of e on the mediator width is highlighted in figure 6 (right panel) 
where we plot the total production cross section times efficiency as a function of m x . Curves 
corresponding to the parameter settings g^ = 1/3 and g^ = 3 (g£ = 1/6 and g x = 3/2) are 
marked in blue (orange). Dashed lines are obtained by assuming the (nominal) efficiency e 
provided in [20], whereas solid lines are calculated using the (actual) efficiency that takes 
into account the broad-width effect on the dijet invariant-mass cut. While for the narrow- 
width scenario (cf. orange curves) the discrepancy between nominal and actual efficiency is 
negligible, for the broad scenario (cf. blue curves) the discrepancy is relevant. For instance, 
for the broad case considered in the figure the ratio between nominal and actual efficiencies, 
r e , increases from around 1.2 to 2 as the invisible decay channel opens up. 

Once r e is determined, it is straightforward to recast the UA2 bound in the model we 
are analysing. The UA2 collaboration explicitly applies its constraint to a sequential Stan¬ 
dard Model (SSM) Z' and presents a 90% C.L. bound on a(pp —> ^ssm) x BR(ZggM —> jj) 
(assuming the efficiency e ~ 0.7) [20]. A parameter point in our model is excluded if the 
predicted total cross section o(pp —> R —> jj) violates the condition 

r e a(pp —> R—> jj) <l 

a(pp —r Z SSM ) x BR(Z S g M jj) 

The region of the Mji-rn x plane excluded by this bound is shown in figure 7 (blue regions 
with solid borders), where different coupling scenarios are assumed in the various panels. 11 
For large couplings g^ the UA2 constraint rules out the mass range 130 GeV < Mr < 
260 GeV. For g ^ ~ 1/4 the weakening in sensitivity of the UA2 experiment at rrijj ~ 
200 GeV leads to an allowed horizontal band. The UA2 constraint is not sensitive to a 
resonance with g// < 1/6 for any value of m x . Moreover, for sufficiently large g x , the UA2 
bound cannot constrain the region with m x < Mr/2. 

5.2 ATLAS lepton-associated dijet analysis 

The ATLAS collaboration has analysed a 20.3 fb _1 data sample searching for dijet reso¬ 
nances produced in association with a gauge boson via pp —>• RZ —>• jj£ + £~ and pp —> 
RW ± ->■ jjl±v processes at yfs = 8 TeV [24] . The presence of at least one lepton in the 
final state allows to highly reduce the QCD background and therefore to probe the region 
with dijet invariant mass 130 GeV < mjj < 300 GeV which is otherwise inaccessible at the 
LHC. We implement these two ATLAS analyses to complement the UA2 constraint in the 
low Mr region. 

We produce pp —>• RZ —»• jj£ + l~ and pp —> RW± —>• jjt^v Monte Carlo events in 
accordance with ATLAS detector specifications. We also impose the cuts enlisted in the 
two analyses. The search for the Z-(!F-)associated production requires p^ ^' 1 > 50 GeV 
where p^^ is the transverse momentum of the dilepton (lepton and $r) system. In 

lx We find the bounds from UA2 to be somewhat stronger than the ones shown in [70], in agreement with 
the analyses presented in [45, 71]. 
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addition, both analyses require two jets with pt > 30 GeV and relative pseudorapidity 
|Ar/yy| < 1.75 (see [24] for details on the cuts). To reconstruct the jets we employ an 
anti-fej [80] algorithm with R = 0.4. As in the UA2 analysis, the final outcome of the event 
analysis is a distribution of dijet invariant masses across several nrijj bins. 

To validate the simulations we produce the dijet invariant-mass distribution for the 
corresponding SM background (mainly Z+jets and VF+jets production) and we find agree¬ 
ment with the one reported in [24]. As a last step, we compare the rrijj distributions 
observed by ATLAS with those predicted by the SM with and without a mediator dijet 
signal. 12 We quantify the (in)compatibility of prediction and observation by means of a 
CL S statistical analysis considering ATLAS measurements on W- and Z-associated pro¬ 
ductions at once. 13 More precisely we determine the confidence level that the signal events 
distributed in all rrijj bins are incompatible with the observed distribution. This statistical 
approach enables us to properly analyse even broad resonances. 

The final result is presented in figure 7 (green regions with dashed borders). We see 
that the ATLAS 95% C.L. bound turns out to be slightly stronger than the one by UA2. 
Notably, the region at Mr ~ 200 GeV where UA2 loses sensitivity for = 1/4 is now 
clearly ruled out. 

5.3 CDF dijet analysis 

The CDF collaboration performed searches for a dijet resonance in pp collision data with 
an integrated luminosity of 1.13fb _1 at yfs = 1.96TeV [21]. This analysis requires two jets 
within the central region of the detector, corresponding to a rapidity \y\ < 1. These jets are 
then reconstructed using a cone jet-clustering algorithm with R = 0.7. While the search 
in principle constrains dijet resonances in the mass range 260 GeV < nijj < 1400 GeV, 
we focus on the region 300 GeV < rrijj < 900 GeV, where the most stringent constraints 
are expected [70]. The lower end of this mass region is the most interesting, because the 
mediator is already too heavy to be produced in sufficient number for the UA2 and ATLAS 
searches above, but is still too light to be disentangled from the SM background in LHC 
dijet searches. 

For our analysis, we mimic the CDF detector characteristics and recast the cuts em¬ 
ployed in the experimental analysis. We validate the analysis by reproducing the SM 
background distribution and checking for agreement with the background reported in [21]. 
To determine the compatibility of the model with the CDF data, we apply an analysis 
analogous to the one discussed above for the lepton-associated dijet search. Given the 
large amount of events at small invariant masses and the various error sources affecting the 
correct description of the QCD background, we also include systematic uncertainties in the 
normalisation of the background in our statistical analysis. The systematic uncertainties 

12 The numbers of observed events in each bin have been digitised from the plot in [24] in order to properly 
estimate the bounds on our theory. Indeed, the bound cannot be directly read off from the exclusion limits 
provided in [24] because the analysis is devoted to a different model. 

13 For the basics of the implementation of the CL S method, see e.g. appendix C of [81] . For the combination 
procedure we follow [82]. 
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Figure 7. Parameter regions excluded by dijet searches at UA2, the Tevatron and the LHC. For 
the left (right) column, we have fixed g = (g^ gf) 1 ^ 2 = 1 (g = 0.5), while the different rows show 
different coupling ratios g^/g^. The UA2 bounds correspond to 90% C.L., all other bounds are 
shown at 95% C.L. 
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have been (globally) fixed in such a way that the CDF bounds can be reproduced when 
considering a narrow resonance. 

The results of our simulations are presented in figure 7. The CDF 95% C.L. bounds 
(orange regions with dotted borders) strongly constrain the MR-m x parameter space. The 
general trend is that for 0.25 < < 1 and 300 GeV < Mr < 700 GeV only the parameter 

space with a large mediator invisible width evades the CDF bound because of a twofold 
effect: for m x < Mr/2 and large enough g x the mediator branching ratio to quarks is 
reduced and its total width broadens, resulting in a loss of sensitivity. 

5.4 ATLAS and CMS dijet analyses 

The CMS and ATLAS collaborations look for dijet resonances in 19.7 fb -1 and 20.3 fb~ 4 of 
data of pp collisions at y/s = 8TeV, respectively [22, 23]. The CMS analysis is somewhat 
more involved than the ones discussed above. It first reconstructs jets by means of an 
anti -kr algorithm with R < 0.5 and selects jets with px > 30 GeV and p < 2.5. Second, 
geometrically close jets (A R < 1.1) are combined into wide jets, which are then selected 
to form a final (wide) dijet system fulfilling \&.pjj\ < 1.5, rrijj > 890 GeV and Hr = 
YhiPr > 650 GeV. The ATLAS analysis is comparably straightforward. It uses the anti -hr 
algorithm with R < 0.5 to reconstruct jets and rejects those with px < 50 GeV. Jets 
are then joined into a dijet system whose leading jets must have rapidity |y| < 2.8 and 
|y*| = 12/iead—2/sublead|/ 2 < 0.6. Finally, only dijets with an invariant mass of rrijj > 250GeV 
are considered. 

We analyse the ATLAS and CMS searches for mediator masses Mr > 500 GeV (for 
smaller masses the LHC loses sensitivity compared to the Tevatron). To validate the 
simulation and data analysis we essentially proceed as explained for the CDF and ATLAS 
lepton-associated searches above. In particular we employ the MCLimit code for both 
ATLAS and CMS results at once, assuming that they are not correlated. We find that 
the ATLAS and CMS 95% C.L. exclusion limits are very powerful (cf. red regions inside 
dot-dashed borders in figure 7) although they can be avoided in the parameter region with 
large invisible width. We also observe from figure 7 that even in the absence of invisible 
decays (i.e. for Mr < 2 m x ) the bounds on Mr become stronger as g ^ is increased. The 
reason is that the enhancement of the R production cross section overcompensates the 
reduction of the detection efficiency due to the broadening of the resonance. 

5.5 Comparison with top-quark pair production at the Tevatron 

Given that searches for dijet resonances yield relatively weak constraints for 300 GeV < 
Mr < 500 GeV, it is interesting to consider other kinds of searches that can potentially 
probe this particular parameter region. Indeed, if the vector mediator has non-zero cou¬ 
plings to the third generation of quarks, its presence can lead to a significant enhancement 
of the total tt production cross section. This cross section has been precisely determined 
by combining CDF and DO measurements at yfs = 1.96TeV with C = 8.8fb -1 [83], leading 
to constraints that are most relevant if Mr is close to the tt threshold. 14 

14 Note that recent LHC measurements of the total tt production cross section, such as [84], do not improve 
on the constraints from the Tevatron, since tt production at the LHC is dominated by gluons. 
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Figure 8. Comparison of the combined bounds from searches for dijet resonances (blue) and the 
constraints from the Tevatron measurements of the total tt production cross section (purple) for 
two particular coupling combinations. 


To obtain the constraints from the Tevatron we simulate tt production at leading order, 
taking into account the interference between the SM contribution and new physics effects. 
If the contribution of the new vector mediator is negligible, we obtain a(pp —> tt) = 6.03pb. 
Comparing this value to recent NNLO+NNLL calculations of tt production at the Tevatron, 
which quote a(pp — > tt) = 7.351 q 33 pb [85], allows us to determine the relevant A'-factor 
to be K = 1.22 ± 0.05. We can then approximately capture the effect of higher-order 
corrections by rescaling all simulated cross sections by this factor. 

The combined results of several tt measurements at the Tevatron determine the tt 
production cross section to be cr (pp —>• tt) = 7.60 ± 0.41 pb [83]. Combining experimental 
and theoretical uncertainties, we can therefore place an upper bound on the rescaled cross 
section of a (pp —> tt) < 8.6 pb at 95% C.L. 

As expected, the resulting constraints are most relevant for mediator masses just above 
the tt threshold, i.e. Mr ~ 400 GeV. Moreover, just like the constraints from dijet 
resonance searches, we find the strongest constraints if either m x > Mr/2 or g ^ 
such that the mediator decays dominantly into quarks. In other words, measurements 
of the tt cross section constrain the same parameter space also probed by searches for 
dijet resonances at the Tevatron and we find that the latter constraint is always slightly 
stronger. Two particular examples are shown in figure 8, where we compare the combined 
dijet bounds from figure 7 (blue) with the parameter region excluded by tt measurements 
(purple). For g/p < 1/6 we do not obtain any constraints from tt measurements. Since the 
constraints from tt do not improve upon the constraints discussed above, we will not show 
them in the remaining plots. 

Nevertheless, it is worth emphasising that bounds on the total tt cross section do 
provide an important independent constraint for large quark couplings. We find that for 
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g‘q = 0.5 (and either small g x or large m x ) the Tevatron can essentially exclude the mass 
range 330 GeV < Mr < 700 GeV, while for g£ = 1 the entire range Mr < 900 GeV is 
excluded. In the latter case, the tt cross section can be as large as four times the SM 
prediction (for Mr ~ 400 GeV), making it clear that this case is solidly ruled out. 

6 Combined constraints 

We are now in the position to combine the various constraints discussed in sections 3-5. 
Our central results are shown in figures 9-11. In this section, we discuss these figures in 
detail and draw our conclusions on the complementarity of the different searches. 

6.1 Universal couplings 

In figure 9 we consider the case that the mediator couples to all quarks with equal coupling 
g the different panels corresponding to different coupling strengths. For the left (right) 
column, we have fixed g = (g x g^) 1 ^ 2 = 1 (g = 0.5), while the different rows show different 
coupling ratios g x /g J q as discussed in section 2. We make the following observations: 

• For <7 = 1 , monojet searches are more strongly constraining than direct detection 
experiments for the entire parameter region where the DSP is underproduced in 
the early Universe (and the theory is perturbative). Considering coupling ratios 
g x /g^ > 1 further suppresses direct detection relative to monojet constraints, be¬ 
cause the larger direct annihilation into pairs of mediators and the smaller mediator 
width reduces the DSP abundance. 

• For g = 0.5, monojet searches have very limited sensitivity for the parameter re¬ 
gion m x > Mr/2, corresponding to off-shell production of the DSP. Direct detection 
experiments therefore typically have an advantage over monojet searches in this re¬ 
gion. Nevertheless, monojet searches are still clearly more constraining than direct 
detection for m x < Mr/2. 

• For g x /gq = 1 and the coupling strengths under consideration essentially the entire 
mediator mass range 130 GeV < Mr < (2.5-3) TeV is excluded by the dijet searches 
discussed in section 5. These constraints can be weakened by considering g x /gf > 1, 
both because the production cross section of the mediator is reduced and its invisible 
branching fraction is increased. In particular, for g'^/g^ = 9 dijet searches are insen¬ 
sitive to the parameter region m x < Mr/2, where the invisible branching fraction of 
the mediator is larger than 80%. In this case, there is a very strong complementarity 
between monojet searches and dijet searches. 

• For g^ > 1/2 the combination of the constraints from above excludes essentially the 
whole range of DSP and mediator masses apart from two special regions, namely 
m x ~ Mr/2 > 1 TeV and rn x S> Mr ~ Mz- For smaller quark couplings and large 
ratios g x /gf additional allowed parameter regions open up for m x > Mr. 
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Figure 9. Combined constraints (at 95% C.L.) from direct detection (orange, dotted), searches for 
monojets (green, dashed) and dijets (blue, dot-dashed) compared to the parameter region excluded 
by DSP overproduction (red) and perturbativity (grey). For the left (right) column, we have fixed 
g = (g£ g'q) 1 / 2 = 1 (g = 0.5), while the different rows show different coupling ratios g^/gf- 
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Figure 10. Same as figure 9 but for alternative coupling structures. The left (right) plot shows 
the case of no heavy-quark couplings (isovector couplings). 

For even smaller couplings than the ones considered in figure 9, the bounds from the 
LHC can be significantly relaxed and we find that there is at present no sensitivity for 
g^, g x < 1/6. At the same time, the constraints from perturbativity are weakened, opening 
up additional parameter space for m x Mr. Nevertheless, direct detection constraints 
remain strong even for such small couplings and it becomes increasingly difficult to avoid 
an overproduction of the DSP. 

6.2 Alternative coupling structures 

In figure 10, we consider two alternative scenarios for the coupling structure between the 
mediator and SM particles, namely couplings only to light quarks ( g^ = g// = 0) and 
isovector couplings (g// = —gf etc.). In both cases, we have fixed g// = g x = 1. 

The absence of couplings to heavy quarks changes the results compared to the case of 
universal couplings (see top-left panel of figure 9) due to three separate effects: 

1. The most significant difference comes from the absence of the annihilation channel 
XX > tt. which becomes important for rn x > mt ■ As a result, the DSP abundance 
will be much larger than in the case of universal couplings and hence relic density 
constraints and direct detection constraints will become stronger. 

2. As a somewhat more subtle effect, the reduced number of decay channels means that 
the mediator will have a more narrow width (for fixed coupling strength), leading to a 
smaller relic density and hence weaker constraints from the relic density requirement 
and direct detection in the resonance region m x ~ Mr/2. 

3. For the same reason, the mediator will have a larger invisible branching fraction, 
enhancing the constraints from monojet searches in the parameter region rn x < Mr/2 
(see also figure 5). 
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Figure 11. Same as figure 9 but in the presence of vector couplings. The left (right) plot shows 
the case where g v g A ( g A 3> g l ). Note the change of scale in the left figure. 


We find that as a result of all of these effects direct detection is enhanced relative to LHC 
searches, although the combined bound is still dominated by the LHC for g = 1. 

The sign difference for isovector couplings will not change LHC cross sections, decay 
widths and relic density calculations, all of which are independent of interference effects. 15 
However, interference between the individual quark contributions is crucial for direct de¬ 
tection. One can immediately infer from the values for A q( N ) given in eq. (3.7) that there is 
destructive interference between the up-quark and the down-quark (as well as the strange- 
quark) contribution if all couplings have equal signs. For isovector couplings, on the other 
hand, the interference is constructive. As a result, this simple modification of the standard 
scenario increases the effective DSP-nucleon coupling by a factor of 3 and hence boosts 
direct detection cross sections by an order of magnitude. In spite of this enhancement, 
however, we still find LHC searches to more constraining than direct detection even for 
isovector couplings. 

6.3 Vector couplings 

To conclude this section, let us briefly consider non-zero vector couplings g 1 and g x . Fig¬ 
ure 11 (left) shows the combined bounds for g A = g£ = 0 and g 1 = g x = 1. As expected, 
direct detection receives a huge enhancement due to the presence of SI interactions. As 
a result, LHC searches are essentially irrelevant in comparison. The only exception is the 
parameter region m x < 10 GeV (not shown in figure 11), where direct detection loses 
sensitivity. This region, however, is already fully excluded by the combination of monojet 
searches and the relic density requirement. 

15 One notable exception are DM searches in the mono-W channel, where the relative sign between up- 
quark and down-quark coupling plays an important role [86]. For isovector couplings, these searches are 
expected to give stronger bounds than conventional monojet searches (see e.g. [87]), but a detailed analysis 
in the context of our model is beyond the scope of the present work. 
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As a somewhat more interesting case, let us consider the case gf,g^^> g\ , g^ > 0- 
For such a setup, the vector couplings give no relevant contribution to the total width of the 
mediator and they play no role for the LHC phenomenology of the model. Nevertheless, 
small vector couplings can potentially give a relevant contribution to the relic density 
calculation (due to the absence of a helicity suppression in the s-wave contribution) and 
they will certainly lead to an enhancement of event rates in direct detection experiments. 
One particular example is shown in figure 11 (right). For g J ^ = g^ = 1 and g 1 = g\ =0.1 
direct detection experiments clearly give stronger bounds than monojet searches, but they 
are less constraining for heavy mediator masses than dijet searches. 

7 Discussion 

As we have seen in the previous section, the parameter space under consideration can be 
tightly constrained by combining results from a variety of colliders searches. This obser¬ 
vation immediately leads to two important questions: First, can the remaining parameter 
space be constrained even further using near-future experiments? And second, can these 
constraints be evaded by modifying the assumptions on the coupling structure and the 
cosmology of the dark sector? We will now address both of these questions in turn. 

7.1 Future prospects 

Within the next decade we can expect rapidly increasing target masses and exposures in 
direct detection experiments as well as a significant improvement in the centre-of-mass 
energy and luminosity at the LHC. To discuss the expected impact of these experimental 
developments, we consider two different projections: 

1. A three-year projection considering the LHC with an integrated luminosity of lOOftV 1 
at 14 TeV and a direct detection experiment based on liquid xenon with a total 
exposure of 2 ton-years (like XENON1T [88]). 

2. A ten-year projection considering the LHC with an integrated luminosity of 300 ftV 1 
at 14 TeV and a direct detection experiment based on liquid xenon with a total 
exposure of 10 ton-years (like LZ [89]). 

For the projected sensitivity of direct detection experiments we assume that these 
experiments will remain essentially free of backgrounds so that the resulting bounds on 
the DSP-nucleon scattering cross section is inversely proportional to the total exposure. 
XENON1T can hence improve the current bound from LUX by a factor of 70, while LZ 
can achieve a factor of 350 improvement. 16 

In [91], the ATLAS collaboration has studied the sensitivity of the LHC at 14 TeV for 
monojets. We have generated sets of monojet events based on the cuts suggested there, 
imposing in particular a very stringent cut on missing transverse energy of > 800 GeV. 

16 As long as we consider only DSP masses above 10 GeV, the sensitivity of LZ is not limited by the 
background from solar neutrinos [90]. 
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For these cuts, the expected bounds on the fiducial monojet cross section are around 1 fb. 1 ' 


It is important to note that while a huge improvement in sensitivity will be achieved from 
the larger centre-of-mass energy, which better separates signal from background, the gain 
from growing luminosity is rather slow due to significant systematic uncertainties in the 
background estimation (assumed to be 5%). 

Various previous works have discussed the projected sensitivity of dijet searches at 
future colliders. Ref. [71], for example, has considered prospects for searches for dijet 
resonances in association with SM gauge bosons. The best expected sensitivity was found 
for the case where the gauge boson is a photon. We take the projected bounds for the one- 
Higgs doublet case (called BP1 in [71]), having confirmed explicitly that the production 
cross section of the mediator is the same for vector and axial couplings. A dedicated study 
of dijet resonances at higher energies has been performed in [92], We adopt their bound for 
an integrated luminosity of 300 fb -1 and infer the corresponding bound for C = 100 ftV 1 * by 
using that bounds on the coupling g q scale with luminosity approximately proportional to 
T 1 / 4 , because bounds on the production cross section should become stronger proportional 
to \JH and the production cross section is proportional to g q . ls 

We note that (projected) constraints on dijet resonances are typically shown as con¬ 
straints on the mediator-quark coupling, i.e. g q < g max (Mn), assuming that the mediator 
can only decay into quarks. In the presence of an additional invisible decay channel, these 
constraints can be considerably weakened. To apply constraints from the literature to our 
model, we therefore have to rescale any bound according to 19 



(7.1) 


Presently there are no sensitivity studies for the intermediate mass region 600 GeV < 
Mr < 1200 GeV and, in light of the large resources required to deal with the huge statistics 
in this mass range in dijet searches, it is not clear whether future colliders with higher 
centre-of-mass energy can significantly improve their sensitivity. An interesting option [93] 
to overcome this difficulty consists of recording either a fraction of the events (this procedure 
has been already applied in [23]) or a reduced part of each event (see e.g. [94]). A different 
possibility to make progress may be searches for tt resonances along the lines of [95, 96], 
which would also be a generic signature of the model considered here. We leave a detailed 
study of these searches for future work and simply use existing bounds from the Tevatron 
and the LHC for the intermediate mass region. 

In figure 12 we consider the case g q = 1/4, g£ = 1 5 which is comparably weakly 
constrained by current experiments (see centre-right panel of figure 9). We can see direct 

1 ' The precise numbers can be worked out by comparing the numbers of predicted events given in table 3 

with the expected bounds shown in figure 6 of [91]. 

18 Ref. [92] considers only the case of very narrow resonances. Nevertheless, we have confirmed explicitly 
that one obtains comparable projections for Tr/Mr « 0.03-0.05, corresponding to the couplings considered 
in figure 12. 

19 An additional rescaling factor of 1/6 arises when comparing with [70, 92] due to different conventions 
for the couplings. 
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Figure 12. Same as figure 9 but showing projected sensitivities for near-future experiments. We 
focus on the case g// = 1/4 and g£ = 1 (cf. centre-right panel of figure 9) and show a three-year 
(ten-year) projection on the left (right). See text for details. 

detection experiments rapidly gaining sensitivity in the parameter region m x > Mr/ 2, 
while monojet searches can probe the resonance region m x ~ uir/ 2 up to mediator masses 
of around 2 TeV. Finally, dijet searches will be sensitive to values of Mr up to (3-4) TeV. 
We note that a possible future 100 TeV hadron collider could easily reach up to Mr ~ 
10-20 TeV, placing extremely strong limits on candidate dark sector theories [34, 92, 97]. 

7.2 Generalised bounds 

Let us now turn to a discussion of the model assumptions that we have made. Crucially, 
all of our results have been based on the assumption that the width of the mediator is 
dominated by the couplings between the mediator and quarks as well as DSPs. Similarly, 
we have assumed that only these couplings are relevant for the relic density calculation. 
While it is justified to neglect additional couplings to leptons and SM gauge bosons given 
the stringent experimental bounds, it is of course conceivable that there are additional 
unstable states in the dark sector, which are light enough to provide additional channels 
for DSP annihilation and/or mediator decay. If these additional light states subsequently 
decay into high-multiplicity SM states, it is conceivable that they could evade detection in 
all existing collider searches. The presence of such additional states implies that we can no 
longer calculate Tr and Ddsp in terms of the couplings g q and g x and have to treat them 
as additional free parameters. 20 

In figure 13 we show two examples for how our bounds would change if we fixed Ddsp = 
Ddm rather than calculating the relic density in terms of g q and g x . By construction, there 
is no longer an excluded parameter region corresponding to the overproduction of the 
DSP, so the parameter region Mr 1 TeV 3> m x is no longer excluded. At the same 

20 In principle, the presence of additional light states can only lead to an increase in IV and a decrease 
in Pdsp, so that these parameters cannot be chosen arbitrarily for given g q and g x . However, we neglect 
this complication for the purpose of the present discussion. 
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Figure 13. Same as figure 9 but fixing the DSP relic abundance Ddsp = ^dm- As a result, there 
is no longer an excluded parameter region corresponding to DSP overproduction and the constraint 
from direct detection is significantly strengthened in the parameter region where previously the 
DSP was underproduced. To make this change of perspective explicit, we show the direct detection 
bound in purple rather than orange (cf. figure 4). 


time, bounds from direct detection are significantly strengthened in the parameter region 
where previously the DSP was underproduced. Consequently, these bounds are now more 
competitive with monojet and dijet searches, which are not affected by changing the value 
of Ddsp- 

In contrast, increasing the total width of the mediator will strongly affect dijet searches, 
because a broader resonance will be less visible in the dijet invariant-mass distribution and 
furthermore the presence of additional decay channels will reduce the branching ratio of 
the mediator into light quarks. The modified branching ratios also imply that there will 
be weaker bounds from monojet searches in the parameter region where the mediator can 
decay into DSPs. For heavy DSP masses, on the other hand, the monojet cross section is 
largely independent of the mediator width (cf. section 4). 

To study the dependence of the monojet and dijet bounds on the mediator width in 
more detail, we show in figure 14 the experimental bounds on x BR(R —> inv) (left 
panel) and x BR(R —> jj) (right panel) for different values of Ti?. 21 The combination 
gq x BR is chosen because — within the validity of the narrow-width approximation - 
the resulting bound is expected to be independent of the mediator width. Indeed, the 
left panel of figure 14 confirms explicitly that for m x -C Mr the monojet cross section is 
proportional to g^ x BR(i? — > inv) even for rather broad resonances, i.e. the dependence 
on the width of the mediator only enters via the branching ratios. 

In the right panel of figure 14, on the other hand, we can clearly see the loss of 
sensitivity of dijet searches for broad resonances. As the mediator width is increased, the 
bound on g^ x BR(R — > jj) becomes weaker. Note that, in the present work we do not 

21 Note that these bounds apply for axial couplings as well as vector couplings. If both kinds of couplings 
are non-zero, the bounds apply on the combination g 2 = (gf) + (g)') ■ 
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Figure 14. Experimental bounds on g % x BR(A —> inv) (left) and g ^ x BR(7? — > jj ) (right) as a 
function of Mr for different values of Tr. 

consider resonances with a width larger than 50% of its mass. For such resonances it is very 
difficult to distinguish the dijet invariant-mass distribution of the signal from backgrounds. 
A more promising strategy to constrain very broad resonances might therefore be to study 
dijet angular correlations [98], such as the ones considered in the context of constraining 
four-fermion operators (see e.g. [99]). We leave this line of investigation for future work. 

To conclude, we emphasise that the bounds shown in figure 14 can in principle be 
applied to any model containing a spin-1 resonance that can decay into quarks and invisible 
particles (where the latter could also be SM neutrinos). For a given mass of the resonance, 
one simply needs to calculate the total width of the resonance and the respective branching 
ratios and then compare the model prediction with the appropriate bounds. The specific 
model of an axial-vector mediator, which we have chosen to present our results, therefore 
only serves for the purpose of illustration, but does not limit the applicability of our 
analysis. 

8 Conclusions 

In this work, we have discussed the phenomenology of a dark sector particle (DSP) with 
relatively large couplings to quarks. We considered an axial-vector mediator as a suitable 
example to illustrate the compelling interplay between different DM detection techniques, 
as well as the impact of additional theoretical considerations. The relevant constraints arise 
from collider searches for missing energy and resonances, DM direct detection experiments, 
the DSP relic abundance and perturbativity, each of which has a unique foothold in the 
search for DSPs. 

To calculate constraints from direct detection experiments, we calculate the DSP relic 
abundance Ddsp based on the assumed interactions and account for the fact that the local 
DSP density is proportional to Ddsp- This approach suppresses constraints from direct 
detection experiments in the parameter region where the DSP is underproduced. Moreover, 
this approach implies that direct detection bounds are largely invariant under a rescaling 
of couplings (see figure 4), making them less sensitive than LHC searches whenever the 
couplings between the DSP and quarks are large. 
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Two classes of LHC constraints are important in the present context: Searches for DSP 
production in events with large amounts of missing transverse energy in association with 
SM particles (e.g. monojet events) and direct constraints on the mediator from resonance 
searches. While the former have been intensely studied in the literature, the latter have 
received much less attention in the context of DM searches. In fact, it has often been 
assumed that these kinds of searches give only very weak constraints on mediators with 
small mass or large width. 

In the present work we provide a comprehensive analysis of searches for dijet resonances 
in UA2, the Tevatron and the LHC across a wide range of mediator masses. In particular, 
we reinterpret an existing ATLAS analysis for dijet resonances produced in association with 
leptonically decaying gauge bosons in the context of our model and obtain strong bounds 
in the region where the new mediator R is rather light, namely 130 GeV < Mr < 300 GeV. 
Furthermore, we show that, although the experimental searches become less constraining 
for broader resonances, there are still stringent bounds on mediators even with width 
r R ~ Mr/ 2 (see figure 7). These bounds need to be taken into account when interpreting 
DM searches at the LHC in terms of simplified models with an s-channel mediator, because 
they apply to a wide range of models and in many cases complement or even surpass other 
search strategies (see figure 14). 

For example, for the choice of (axial) couplings of R to quarks and DSPs g// = g x = 1 
all mediator masses in the range 130 GeV < Mr < 3 TeV are excluded by dijet searches, 
independent of the DSP mass. Only two small parameter regions remain viable, namely 
the resonance region Mr ~ 2 m x 1 TeV and the low-mass region m x 2> Mr ~ Mz , 
where the DSP annihilates directly into the mediator (see figure 9). These constraints 
can be significantly weakened by considering g x /gf 1 and we have studied the cases 
g x /g'q = 4 and g x /gf = 9 in detail. For example, for g£ = 1/4 and g x = 1 there is a 
strong complementarity between monojet and dijet searches, but there are still significant 
mass ranges unconstrained by present data. These couplings thus provide an interesting 
benchmark point to study projections for near-future experiments (see figure 12). 

We have furthermore pointed out that in the presence of non-zero axial couplings, it 
is not possible within a simple perturbative theory to raise the DSP mass rn x arbitrarily 
above the mediator mass Mr. The same problem also manifests itself in the annihilation 
process XX tt violating perturbative unitarity for large DSP masses. This observation 
implies that a significant part of the MR-m x parameter plane is theoretically inaccessible. 
In combination with the constraints on the mediator mass from the relic density calculation, 
we obtain upper bounds on both the mediator mass and the DSP mass for fixed couplings, 
so that the relevant parameter space is necessarily closed. 

With improving sensitivity at direct detection experiments and colliders it will be pos¬ 
sible to probe the allowed mass ranges for smaller and smaller couplings, up to the point 
where DSP annihilation into quarks will generally be insufficient to avoid DSP overproduc¬ 
tion. In other words, there is the realistic chance to comprehensively test the idea that a 
DSP can have large interactions with SM quarks mediated by a spin-1 particle. Such a con¬ 
clusion will not be achieved by a single experiment, but it will require significant progress 
across a range of different experimental strategies combined with continuing theoretical 
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studies of less-explored search channels. 
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